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Abstract 

Self-propulsion of cellular microswimmers generates flow signatures, commonly classified as 
pusher- and puller-type, which characterize hydrodynamic interactions with other cells or bound¬ 
aries. Using experimentally measured beat patterns, we compute that flagellated alga and sperm 
oscillate between pusher and puller. Beyond a typical distance of 100 /im from the swimmer, inertia 
attenuates oscillatory micro-flows. We show that hydrodynamic interactions between swimmers 
oscillate in time and are of similar magnitude as stochastic swimming fluctuations. 
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For a single cell swimming in a fluid, inertia is negligible mq. Cellular swimmers ex¬ 
ploit propulsion strategies independent of inertia that allow for net propulsion by viscous 
forces [3]. Periodic, non-reciprocal body shape changes I3H7] conquer the time-reversibility 
of the Stokes equation, which governs fluid flow in the inertia-less limit of zero Reynolds 
number. For example, traveling bending waves of long slender cell appendages termed flag¬ 
ella propel many eukaryotic cells, including sperm or swimming alga j8]. Self-propulsion by 
such periodic shape changes implies that net motion is a second-order-effect, superimposed 
to an oscillatory motion, first-order in the amplitude of the swimming stroke pun]. 

The flow fields generated by microswimmers are commonly classified as pusher- or puller- 
type, depending on the direction of fluid flow along the axis of net swimming [TT1 - IT3] : 
inward flow towards the swimmer characterizes a puller, while outward flow characterizes 
a pusher, see Fig. [^b). This characterization allows a simple assessment of the interaction 
of microswimmers with boundary surfaces or inter-swimmer interactions. For example, 
it was suggested that pushers become hydrodynamically attracted to boundary surfaces 
by an inward flow perpendicular to the axis of net swimming [T]. Swimmer type further 
determines active rheological responses such as shear thinning in dense suspensions of micro¬ 
swimmers [nils]. 

Here, we compute time-varying flow fields induced by flagellated swimmers, using exper¬ 
imentally measured beat patterns [iniiiz]. We show that flagellated microswimmers peri¬ 
odically oscillate between pusher- and puller-type swimming, which implies more complex, 
dynamic interactions with boundaries and other swimmers. Further, we discuss how inertial 
effects, usually neglected in microswimming problems, attenuate oscillatory micro-flows at 
a distance S ~ 100 pm from the swimmer. For that aim, we present a novel approxima¬ 
tion scheme to account for inertial effects in microswimming. From our analysis, we And 
that not only do flow fields oscillate, but also the rate of hydrodynamic dissipation associ¬ 
ated with flagellar swimming. This simple fact has important consequences [IE]. It implies 
that flagellar beat patterns are not optimized to minimize hydrodynamic dissipation as sole 
optimization criterion, as considered previously [T91 - I2T] . 

a. Time-dependent micro-flows. Swimming by periodic body shape changes results in 
time-varying flow fields, which have been measured experimentally for a flagellated model 
swimmer, the bi-flagellated green alga Chlamydomonas [22] • Using experimentally mea¬ 
sured beat patterns [H], we computed time-dependent flow fields, first for the limit of 


2 


zero Reynolds number, see Figure [^a,b). Specifically, we used a fast boundary element 
method [2H] to determine the surface density of forces f(r, t) exerted by the time-dependent 
surface S(t) of the swimmer. The dynamics of S{t) is a superposition of an imposed shape 
change of the two flagella, and a rigid body motion of the whole swimmer, which is solved 
self-consistently to ensure force and torque balance of the whole cell. 

In the limit of zero Reynolds number, the fluid flow field v(r, t) is obtained by propagating 
f(r, t) with the Green’s function G of viscous flow, Uj(r, t) = Gij{r—r')fj{r', t). Here, 

Gij{r) = {5ij + rjrj)/(87rr]r), is known as the Oseen tensor, rj is fluid viscosity, and r = |r|, 
r = r/r. Next, we show how this flow field can be decomposed into fundamental singularities 
with an oscillating force dipole characterizing the far field. 

b. Hydrodynamic multipole expansion. The force density fj{r) can be decomposed into 
Cartesian multipoles (r — TqY fj{r,t) [211, similar to the decomposition of 

a charge distribution into charge multipoles. Here, J = (R,... ,jk) denotes a multi-index 
and we use standard multi-index notation. The force density is formally recovered from its 
multipoles as fj{r) = ^j(—— ro)/J! [2H]- We can thus represent the flow 
field as a superposition of fundamental singularities 

«’.('•) = ( 1 ) 

These singularities are traditionally known as Stokeslet (| J| = 1), Stokes doublet (|J| = 
2) m- Note that a self-propelled swimmer, being free from external forces and torques, 
does not exert any net force or torque on the surrounding fluid. Thus, the force monopole 
Hj vanishes, and the force dipole matrix is symmetric. This force dipole represents 
the leading order singularity, which characterizes the far field of fluid flow, decaying as r~^ 
with distance r from the swimmer. For an incompressible fluid, the trace of does not 
contribute to the flow field, but describes a localized pressure. Without loss of generality, 
we assume that Hj^k is traceless. 

The force multipoles change in time, as the microswimmer changes shape and position. 
The leading order singularity Hj^k oscillates with the beat frequency oj [37] . 

We find that higher order multipoles dominate the near field on length-scales comparable 
to the size of the swimmer, while the dipole contribution faithfully reproduces the far field 
at distances of 10 — 100 pm, see Fig. [^a-c). At even larger distances, inertial effects cannot 
be neglected anymore, as we discuss next. 
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FIG. 1: Time-dependent fluid flows induced by a Chlamydomonas cell swimming in water. (a,b) 
Near field and far field of fluid flow, computed in the limit of zero Reynolds number, corresponding 
to the idealization of an inertia-less fluid, (c) The far field is faithfully reproduced by its leading 
order singularity, the Stokes doublet flow induced by a localized force dipole, cf. Eq. 0. (d,e) 
Beyond a characteristic distance 5 from the swimmer, inertial effects attenuate oscillatory micro¬ 
flows, resulting in a deviation from the Stokes doublet flow. Shown are solutions of the linearized 
Navier-Stokes equation, Eq. Q, induced by the oscillating force dipole from panel (c). The white 
circle in (e) has radius 35. Beat frequency: lo/{ 2Tr) — 50 Hz. 

c. Screening of oseillatory miero-flows by fluid inertia. The full nonlinear Navier- 
Stokes equation of fluid flow contains both viscous and inertial terms pv -|- pv • Vv = 
—Vp-l-pAv. The inertial terms on the left-hand-side are usually negligible in the proximity 
of a microswimmer. Their relative magnitude compared to viscous terms is approximated by 
the instationary Reynolds number, Re^^ = pu)A?/r}. For a typical flagellated swimmer in wa¬ 
ter, Re^, is small: Re^, = 10“^ for amplitude A = b pm, beat frequency a;/(27r) = 50 s“^, fluid 
density p = lO^kg/m^, viscosity p = 10“^ Pa-s. Yet, inertial effects become important at a 
distance. For the flow field induced by a force dipole, the unsteady acceleration pv ~ pu)r~^ 
decays slower than the viscous forces pAv ~ r“^, while the convective acceleration decays 
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like pv • Vv ~ r Inertial effects due to unsteady acceleration are thus expected to become 
important beyond a distance 5 = [2ri/{pu))]^/^ [SH]- Using the above values, we estimate 
S 100 /im. 

This argument further shows that convective acceleration can be neglected throughout 
the fluid. Hence, the linearized Navier-Stokes equation applies 


pv = — Vp + pAv. 


( 2 ) 


The flow induced by an oscillating force had been studied by Stokes more than a century 
ago [26l [27]. In modern language, an oscillating force monopole fjexp{—iut) induces an 
oscillating flow field Vi{t)(r) = exp(—icnt) characterized by the instationary Stokes 

tensor 

Gij{r,u) = [(1 - iV)5ij + (1 + 3tT>)fifj], (3) 

where £ = exp[—(1 — i)r/S] and V = {8/rY[£~^ — 1 — (1 — i){r/S) — (1 — z)^(r/(5)^/2]. 
For a Fourier transform, see [28|. The near field with r -C 5, reproduces the Oseen tensor 
Uj(r) = Gij{r)fjexp{—iut) + 0{r/5). For the far field with r S> 5, we obtain a potential 
flow that decays as ~ r~^ with Uj = —3i5AGij{r)fjeyip{—iu)t) + 0(5/r)^. The rotational 
part of the Stokeslet flow field is exponentially attenuated as exp(—r/5). In analogy to the 
common name ‘Stokeslet’ for flow fields of the form Gjjfj, one may call the oscillatory flow 
field Qijfj exp{—iu)t) ‘oscilet’. 

In the limit S ^ L, we can estimate the oscillating force dipole exerted by the 
microswimmer using the limit of zero Reynolds number as employed above. The far field 
according to Eq. Q, including inertial effects due to unsteady acceleration, is then obtained 
as a superposition of oscilet derivatives, 


exp(-mn;t). 


(4) 


where exp(—za;nt). Figure |^d,e) show resultant flow fields for Chlamdy- 

omonas, exemplifying the effect of inertial screening. 

d. Oscillating between pusher and puller. We introduce a coordinate system (ei, 62 , 63 ) 
of ortho-normal vectors, slowly co-moving with the microswimmer, such that ei points along 
the net swimming direction. With this choice of coordinates, the component J^n of the force 
dipole tensor determines the flow speed along the swimming direction as vi = 3 J^ii/( 87 rpr^) 
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(in the limit of zero Reynolds number, neglecting higher multipoles). Thus, the force dipole 
tensor allows to discriminate pusher-type (Jhi > 0) and puller-type (Jhi < 0) swimmers. 
More generally, we can write the traceless, symmetric force dipole tensor as 


= fs^ 



^2 0 0 ^ 


^0 0 0 ^ 



0-10 

+ 

0 

0 



1° 0 -ij 


^0 0 —Cy 





(5) 


where fa denotes a scalar dipole strength, e a dimensionless asymmetry parameter, and 
a rotation matrix. 

For a microswimmer with special symmetries, ^)^ can be chosen particularly simple. The 
idealized bacterium shown in figure [^a) has rotational symmetry around the ei-axis, and 
we may choose = I and e = 0. We find /* > 0 , confirming this swimmer as a pusher. 

Next, for a swimmer with double mirror-symmetry, such as Chlamydomonas with sym¬ 
metric breast-stroke beat, we can choose Tl=][. In this case, the sign of fa determines 
pusher-type (/« > 0) or puller-type (/« < 0). Remarkably, we find that the sign of fa alter¬ 
nates during a beat cycle, being negative during the effective stroke, but positive during the 
recovery stroke, see Fig. [^b). Further, the asymmetry parameter e is found to transiently 
exceed unity, implying that the flow in the ( 02 , e 3 )-plane normal to the swimming direction 
changes sign as a function of azimuthal angle. We will use the terms ‘asymmetric pusher’ 
and ‘asymmetric puller’ for a microswimmer with |e| > 1. Note that a smooth transition 
from a pusher to a puller generally passes through an asymmetric swimmer-type if 
The time-averaged force dipole exerted by Chlamydomonas corresponds to an asymmetric 
puller. 

Finally, for a swimmer with single mirror-symmetry, say with respect to a plane with 
normal 63 , we can choose as a rotation around the 63 axis by an angle a. This case 
applies to sperm cells with planar beat patterns ca. swimming in the (ei, e 2 )-plane. Now, 
d'li = /s[2 cos^ cr -|- (e — 1) sin^ «]. Figure |^c) shows the time-dependent force dipole exerted 
by a sperm cell, revealing again an oscillation between pusher- and puller-type swimming. 

e. Hydrodynamic interactions. Microswimmers can interact at a distance with other 
swimmers or a boundary wall by virtue of the flow they generate [29j . It has been proposed 
that hydrodynamic interactions account for the mutual alignment of several swimmers, beat 
synchronization in collections of flagella [30], or the active accumulation at boundary sur¬ 
faces |SI1IS2|. The oscillatory nature of microflows generated by flagellated swimmers implies 
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FIG. 2: (a) A bacterium propelled by a rotating helical filament represents a symmetric pusher, 

(b) Chlamydomonas oscillates between puller- and pusher-type swimming. From top to bottom 
(fifth column shows time-averages): Flagellar shapes at equidistant phases of the beat cycle. Radial 
projections of the induced Stokes doublet flow for each of these phases, cf. Eq. Q. Instantaneous 
swimming velocity as a function of flagellar phase. Dipole strength cf. Eq. i)- Flow signatures 
are classified as either symmetric or asymmetric pusher- or puller-type swimmers, depending on 
whether the rescaled asymmetry parameter (1 — |e|)/(l-|-|e|) is positive or negative, (c) A swimming 
sperm oscillates between pusher and puller, with a force dipole whose orientation rotates during 
the beat cycle. 

that these interactions likewise oscillate: Fig. shows periodic variations in the distance be¬ 
tween two Chlamydomonas cells swimming side-by-side. These interactions strongly depend 
on the relative phase difference of their flagellar beat. While oscillatory flows dominate 
instantaneous interactions, their net contribution can be small compared to that induced by 
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FIG. 3: Hydrodynamic interactions between microswimmers are weak and oscillatory. Sidewards 
motion of a Chlamydomonas cell induced by hydrodynamic interactions with a second cell that 
initially swims parallel at a distance d — 24 /im, with swimming stroke either in-phase or anti¬ 
phase, respectively. For comparison, we show the expected r.m.s. side-wards displacement of a 
single cell due to passive thermal diffusion (orange shading), or active amplitude ffuctuations of 
the ffagellar beat [33] (gray). 

the static component of flow. In fact, the interaction between two static dipoles separated 
by a distance d scales as while the net interaction between two oscillating dipoles scales 
as after averaging over an oscillation cycle. This is because the flow field induced by 
an oscillatory Stokes dipole ^) has to couple to the flow gradient i^d ^) to yield a net 
effect. Further, hydrodynamic interactions can be masked by thermal and active noise, see 


Fig-11 


/. Phase-dependent rates of hydrodynamie dissipation. The viscous flow generated by 
a swimming cell causes continuous dissipation of energy. For the flow fleld induced by an 
oscilet of strength /o cos (tut), we find an oscillating rate of hydrodynamic dissipation outside 
a spherical cut-off region of radius Rq centered around the oscilet 




where Ri = (6/5)5. In fact, this value equals the viscous dissipation associated with a 
Stokeslet flow field induced by an oscillating force monopole /o cos (cut) inside a spherical 
shell bounded hy Rq < r < Ri. Similar results are found for higher multipoles. We conclude 
that the limit of zero Reynolds number is appropriate for computing rates of hydrodynamic 
dissipation in the limit L <C 5. 

Thus, neglecting inertial effects, the instantaneous power output of the swimming cell, 
f(r)-v(r), equals the rate of hydrodynamic dissipation in the bulk of the 
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fluid 2ri f Vv(r):Vv(r). Importantly, any swimming stroke that minimizes the mean 
power output (7Z) = f dip'R{(p) must satisfy 7^((y?)=const [TH]. Otherwise, reparameterizing 
(p'{if>) = 2w d61Z{6)~^/‘^ reduces (JZ). 

We computed 7?.((p) for a swimming Chlamydomonas cell, see Fig. We find a pro¬ 
nounced phase-dependence of the rate of hydrodynamic dissipation, corresponding to a 45% 
higher value for {TV) compared to a beat with reparametrized phase. This implies that 
flagellar beat patterns are not optimized for minimal hydrodynamic dissipation. Interest¬ 
ingly, hydrodynamic dissipation is maximal during the flagellar recovery stroke, for which 
the distance between flagella and cell body is small, implying high local shear rates. Phase- 
dependent rates of hydrodynamic dissipation correspond to phase-dependent active driving 
forces |T6] , which have been proposed to facilitate hydrodynamic synchronization in flagellar 
pairs, e.g. in Chlamydomonas [34] . 

Previously, dissipation rates had been estimated from experimentally measured flow fields 
[22] • These measurements represent hydrodynamic dissipation in the far held, as it is inher¬ 
ently difficult to resolve localized shear flows at distances of a few microns from the swimmer. 
We therefore computed hydrodynamic dissipation in the far held, using a cut-off distance of 
6 pm, obtaining good agreement with experiment [22] . 



FIG. 4: Hydrodynamic dissipation rate varies dnring the beat cycle. We show the total rate 
of hydrodynamic dissipation for a swimming Chlamydomonas cell as a function of flagellar phase 
{TZtoti red), and similarly the rate of dissipation in the far-field (T^far, blue), corresponding to its 
oscillating force dipole (with cut-off at 6 pm). Insets show the surface density of work exerted by 
the cell on the fluid, which sums up to T^-tot- 

g. Conclusion. Using experimentally measured beat patterns, we computed time- 
dependent flow fields generated by flagellated microswimmers and found that these oscillate 
between pusher- and puller-type swimming. We presented a novel approximation scheme to 
account for inertial effects in microswimming problems, and And that fluid inertia attenuates 
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oscillatory flows beyond a typical distance of 100 //m from the swimmer. Within this dis¬ 
tance, oscillatory flows dominate any static component. These flow signatures characterize 
hydrodynamic interactions between several swimmers, which are likewise found to oscillate 
in time. The static component of hydrodynamic interactions can be of similar magnitude 
as both thermal and active fluctuations of flagellar swimming. 

Active swimming implies a continuous dissipation of energy into the fluid, which we 
estimate as 300 pN //m or 7.5 • 10^ ksT per beat cycle for a swimming Chlamydomonas cell. 
This energy equals the chemical potential of about 3 • 10^ ATP molecules j35], which serve 
as chemical fuel for the 3 • 10^ molecular motor domains that power the flagellar beat [36] . 
This implies that either only every tenth motors takes a step during each beat cycle, or that 
the efficiency of energy conversion is low. 

Interestingly, we find that the rate of hydrodynamic dissipation varies during the beat 
cycle. This implies that flagellar beat patterns do not minimize hydrodynamic dissipation as 
sole optimization criterion. We speculate that design constraints of flagellar beat generation 
or internal dissipation within the flagellum might also played a role in the evolution of beat 
patterns. 
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